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. Mid-infrared images from the Spitzer Space Telescope Galactic Legacy 

. Infrared MidPlane Survey Extraordinaire program reveal that the infrared 

source IRAS 03063+5735 is a bowshock nebula produced by an early B star, 

7— i . 2MASS 03101044+5747035. We present new optical spectra of this star, classify 

J> . it as a B1.5 V, and determine a probable association with a molecular cloud 

^ , complex at Vlsr=— 38 - —42 km s in the outer Galaxy near £ = 140.59°, 

b = —0.250°. On the basis of spectroscopic parallax, we estimate a distance 
of 4.0+1 kpc to both the bowshock nebula and the molecular complex. One 
plausible scenario is that this a high-velocity runaway star impinging upon a 
molecular cloud. We identify the H II region and stellar cluster associated with 
IRAS 03064+5638 at a projected distance of 64 pc as one plausible birth site. 
The spectrophotometric distance and linkage to a molecular feature provides an- 
other piece of data helping to secure the ill-determined rotation curve in the 
outer Galaxy. As a by-product of spectral typing this star, we present empirical 
spectral diagnostic diagrams suitable for approximate spectral classification of 
O and B stars using He lines in the little-used yellow-red portion of the optical 
spectrum. 
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Introduction 



Short-lived massive stars typically lie within the OB associations in which they were 
born. However, a significant fraction (17%-50%) having high (30-200 km s" 1 ) space veloc- 
ities wander far from such associations and a re known as "runaways" ( )Blaauwlll96ll ; IStone 
19821 : IConti et al.lll977l ; iGies fc Bolton! Il986[ ). While many runaways are detected directly 



on the basis of large radial velocities or high proper motions, others are inferred by the pres- 
ence of a tell-tale bowshock produced b y their supersoni c moti on thr ough the surrounding 



mediu m. Ivan Buren fc McCrayl (119881 ) , Ivan Buren et al.l (119951 ) , and iNoriega-Crespo et al. 



( 119971 ) used Infrared Astronomical Satellite (IRAS) images to locate arc-shaped bowshock 
features associated with high-velocity O stars. Bowshocks associated with high-velocity O 
stars have been detected more easily using the recent generation of mid-infrared surveys 
with higher angular resolution performed by the Midcourse Space Experiment (MSX) and 
the Spitzer Space Telescope (SST). Recent stud i es of b o wshocks around massive s tars in clude 
Brown fc Bomand (l2005r)iComer6n fc Pasqualil (I2007T) . iGvaramadze fc Bomand (120081 ). and 
Povich et al.l ( 120081 ) . In iKobulnicky et al.l (120101 ) we used the SST Cygnus-X Legacy Survey 
data to identify several probable runaway stars and their bowshocks in the vicinity of the 
Cygnus OB2 association while introducing a novel technique to exploit bowshock physics in 
order to estimate mass-loss rates from massive stars. 

In this contribution we report the discovery that infrared source IRAS 03063+5735 is 
a bowshock nebula associated with an early B star in the outer Galaxy near I = 140.59°, 
b = —0.250 °. Data from the SST Galactic Legacy Infrared MidPlane Survey Extraordinaire 
360 (GLIMPSE 360) mid-infrared survey reveals the bowshock morphology for the first 
time. We use optical spectroscopy from the Wyoming Infrared Observatory (WIRO) 2.3 m 
telescope to determine a spectral type and radial velocity. These data, in turn, enable an 
estimate for the distance to the star and an associated molecular cloud complex. 



2. Data 

2.1. Discovery and Infrared Imaging Data 



GLIMPSE 360 is a SST warm mission legacy program survey of the outer G alaxy in the 
[3.6] and [4.5] /im bandpasses (bands 1 & 2) of the Infrared Array Camera (IRAC lFazio et al. 
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2004 ). As part of a systematic study of "intermediate-mass" star-forming regions where stars 
up to ~B0 are being formed, our team has been investigating IRAS sources having colors 
characteristic of soft UV radiation fields from early-to-mid B stars. IRAS 03063+5735 was 
initially selected as having IRAS colors indicating a such a soft radiation field similar to 



Gala ctic star forming regions hosting intermediate-mass stars (jKertonll2002l ; lArvidsson et al. 



20lM Our visual examination of IRAS 03063+5735 in GLIMPSE 360 images revealed a 



striking bilateral symmetric bowshock morphology at [3.6] and [4.5]. Figure [T] (right panel) 
shows a color representat ion of the nebula w ith the SST [3.6], [4.5], and Wide-Field Infrared 
Survey Explorer (WISE; IWright et al.ll2010l ) band W3 [12.1] fj,m image in blue, green, and 
red, respectively. The nebula has a diameter of ~2'. A star, 2MASS 03101044+5747035 
(hereafter Star A), lies near its center and along the axis of symmetry, making it the probable 
progenitor star. The left panel shows the 2MASS (jThe Two Micron All Sky Survey! 120061 ) 
J/H/K s image in blue/green/red, respectively. Note the presence of two infrared-bright 
rims in the mid-IR images. The first, with an apsis located ~7" from Star A, is yellow-green 
in color, indicating relatively strong emission in the [4.5] (green) band. Diffuse emission in 
this band is typically dominated by shocked molecular hydrogen often seen in ou tflows from 
young stellar objects ( Reach et al. 2006 ; Davis et al. 2007 ; Cyganowski et allbood ). although 
the Brackett a line is also in this bandpass. The second, with apsis ~27" from Star A, is 
relatively brighter in the [3.6] and [12.1] bands, giving it a reddish-white hue in Figure (U, but 
it is present at [4.5] as well. The [3.6] band contains broad polycyclic aromatic hydrocarbon 
(PAH) features while the [12.1] band contains both PAH features and a contribution from 
hot dust grains. 



2.2. WIRO Optical Spectroscopy 

We obtained optical longslit spectra of Star A at the WIRO 2.3 m telescope on the nights 
of 2011 October 23 & 27 using the 900 1/mm grating in first order over the wavelength range 
4300 - 7300 A. A total of seven 600-second exposures were obtained with a 2" x 150" slit 
in ~2" seeing with an east-west slit orientation on 23 October and a north-south slit on 
27 October. Spectra were wavelength calibrated using CuAr lamp exposures obtained after 
each exposure of Star A. The wavelength calibration has an rms of 0.05 A and the spectra 
have a resolution of 3.9 A FWHM based on CuAr lines. Individual spectral exposures were 
reduced using standard flat fielding and wavelength rectification procedures before being 
transformed to the Local Standard of Rest (LSR) velocity frame and combined. 



1 As of this writing, the SIMBAD database lists it as a galaxy, presumably on the basis of its infrared 
colors, which include some star-forming galaxies. 
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We also obtained single 5-10 s exposures of 22 bright B-type dwarfs and giants from 
the Smithsonian Astrophysical Observatory Star Catalog for use in calibrating the spectral 
type of Star A. The spectra were reduced and calibrated in the same manner as described 
above. Table [T] lists these stars along with equivalent width measurements for prominent H 
and He lines. These spectra have typical signal-to-noise ratios (S/N) exceeding 200:1 in the 
middle of the spectral range observed. 



3. Analysis 



3.1. Optical Spectra and Spectral Type of Star A 

Figure [2] shows the average spectrum of Star A between 5300 A and 6600 A. The peak 
signal-to-noise ratio near the middle of the spectrum is 65:1, but it is significantly lower 
toward the blue end beyond the pictured range. Labels denote key spectral diagnostic lines 
including He II A5411 (not detected), and the detected lines of He I A5876 and H<x H/3 and 
H7 are also detected but not labeled owing to poorer S/N. All other spectral features are 



inters tellar, including several diffuse interstellar band features (DIB; iJenniskens fc Desert 



19941 ) and the Na I doublet, AA5889,5895. The presence of He I A5876 indicates a relatively 
hot star, while the apparent absenc e of He II A5411 m eans that Star A is not among the 
hottest O stars and is B0. 5 or later ( IKerton et al.lll999l ). 



We used the robust curve fitting package MPFIT ( IMarkwardtl 120091 ) in IDL to measure 
equivalent widths and uncertainties for the He II A5411, He I A5876, and Ha lines in this 
portion of the spectrum. Initially line positions, widths, and depths were allowed to vary as 
free parameters. He II A5411 is not detected, so we measured an upper limit by fixing the 
width of this line to the value computed for the A5876 line and constraining the position to 
lie within 1 Angstrom (55 km s _1 ) of the rest wavelength. Table [T] gives the measured He 
and H line EWs and uncertainties for Star A. No emission lines were detected either at the 
location of the star or at locations along the slit covering the bowshock nebula. 

The red spectral range covered in our data is not one commonly employed for spectral 
typing of hot stars owing to the paucity of spectral features. Consequently, we developed a 
series of diagnostic diagrams using the EW of He and H lines to measure the spectral type 
of Star A more precisely. Figure [3] shows the ratio of He II A5411 equivalent width to that of 
He I A5876 as a function of temperature/spectral type. C rosses, asterisks, and pluses denote 
early type stars from our measurements of spectra in the IJacoby et al.l (119841 ) spectral atlas 
for luminosity classes I, III, V, respectively. The upper ordinate gives the spectral type of 
these stars. Open squares and triangles mark the measurements of Tlusty (ILanz &: Hubeny 
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20031 ; iHubeny &: Land Il995l ) model atmospheres with solar metallicity and temperatures 
as given on the lower ordinate. We adopt the relation b etween spectra l type /luminosity 
class and effective temperature/surface gravity provided by lMartins et al.l (120051 ) in order to 
ensure consistency between the upper and lower ordinates. The atlas O star dwarfs (pluses) 
and model O star atmosphere dwarfs (triangles) show good agreement. A dashed line marks 
a third-order polynomial fit to the pluses and asterisks, and this expression is given near the 
top of the plot, where t 3 is stellar effective temperature in thousands K. Evolved stars and 
models with lower surface gravity characteristic of evolved stars lie systematically above this 
line. 

Star A has a ratio EW(He II A5411)/EW(He I A5876) < 0.14, making it cooler than any 
star on Figure [31 and therefore, later than B0. We constructed a two additional diagnostic 
diagrams appropriate for cooler, B-type stars, using the line measurements from SAO stars 
in Table [IJ supplemented with a few additional early B star spectra from the Cygnus OB2 
radial velocity survey of iKiminki et al.l ( 120071 . 120091 1 . He and H line EWs were measured 
in the same way using the MPFIT code within IDL. Figure H] plots the EW(HeI A5876) 
versus temperature/spectral type (upper panel) and the ratio EW(Ha)/EW(HeI A5876) ver- 



sus temperature/spectral type (lower panel). As in Figure EJ pluse s denote dwarfs 



Lum. 



V), asterisks denote giants (Lum. Ill), and triangles denote Tlusty (lLanz &: Hubenyl 120071 ) 
model atmospheres having log(g)=4.0 appropriate to B dwarfs. Although the dispersion at 
any given spectral type is considerable, the data exhibit a well-defined monotonic trend. 
The solid line in each panel denotes the measurement for Star A, and dotted lines indicate 
the la uncertainties. The dashed curve is a 3rd-order polynomial fit to all of the pluses 
and triangles. Analytic expressions for these curves appear within each panel, where £3 is 
th e temperature in thou sands K. We match temperature to spectral type using the r e lation s 
of ISchmidt-Kalerl (119821 ) which are very similar to the calibration of lUnderhill et al.l ( 119791 ). 
While the measurements on real stars and model atmospheres compare favorably in the up- 
per panel, the models are systematically displaced toward larger EW(Ha)/EW(HeI A5876) 
ratios in the lower panel. We believe that this is the result of broad Balmer wings which are 
measurable in the model Ha lines but not in the lower S/N data. Interestingly, the giants 
and dwarfs show good agreement in both panels. 

Figure H] shows that the spectral type of Star A is near B0. 5 on the basis the EW(HeI 
A5876) measurement (upper panel) and nearer B3 on the basis of the EW(Ha)/EW(HeI 
A5876) measurement (lower panel). Unfortunately, there are few spectral diagnostics sensi- 
tive to temperature for early B stars in the spectral range available. Traditional spectral type 
diagnostics in the blue end of the optical regime are not helpful in our spectra because of 
low S/N — a consequence of reduced instrumental sensitivity and high interstellar reddening 
which suppresses the blue flux from Star A. We adopt a spectral type of B1.5±1.0. H or He 
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emission lines, often seen in B supergiants, are not evident in our spectrum of Star A. Fur- 
thermore, the Ha and H/3 linewidths in star A are 15-18AFWHM- more typical of dwarfs 
and much broader than the Balmer lines in any of our comparison B giants. Therefore, we 
adopt a luminosity of class of V. 



3.2. Star A Spectral Energy Distribution and Distance 



Table |2] lists eight optical and infrared photometric measurements for Star A compiled 
from available data bases. These include the B2, R2, 12 entries from the USNO-B1.0 catalog 
flMonet et al1l2003h . the JHK entries from the 2MASS Catalog, and SST [3.6] and [4.5] data 
from the GLIMPSE Point Source Catalog. Magnitudes have been converted to Jy using 
standard zero points available in the respective documentation. The optical data points are 
particularly uncertain, as they come from photographic surveys, so we adopt uncertainties of 
15%. The infrared photometry from 2MASS and GLIMPSE has uncertainties of 0.03 mag, 
or about 3%. The red colors in any two bandpasses suggest a high degree of reddening for 
an early B star which has close to zero intrinsic color. 

Figure [5] shows the optical and infrared spectral energy distribution of Star A. The 
solid curve is a blackbody approxi mation to a B1.5 V star with temperature T=23,000 K, 
radius R=6.1 (interpolated from Drilling &: LandofifeOOOh . and A v =5.6 mag of reddening 



an 



ac cording to a ICardelli et al.l (119891 ) extinction curve in the optical and near-infrared and 



Indebetouw et al.l (120071 ) reddening law in the mid-infrared. The data are well fit by the 
model for a distance of ~4.0 kpc. The fit is quite good in the infrared where the photometry 
is most accurate and less good in the optical regime where photomteric uncertainties are 
considerable. Nevertheless, the SED constrains the extinction with within about AAy rs 0.4. 
No infrared excess is apparent, as might be expected for evolved massive stars having free- 
free emission or stars hosting circumstellar material. This is consistent with our adoption 
of a dwarf (i.e., main sequence) luminosity class. Based on the fit in Figure we adopt a 
distance of 4. Oil kpc to Star A and the bowshock nebula, where the uncertainty is dominated 
by the uncertainty of the star's spectral type and intrinsic luminosity. At this distance 1" 
corresponds to 0.019 pc and the ~2' angular diameter of the bowshock nebula corresponds 
to 2.3 pc. This also means that the 7" angular distance from Star A to the bowshock 
apsis translates to a linear distance of >0.13 pc, where the inequality accounts for the 
unknown inclination relative to the line of sight. Given that bowshock nebulae are most 
easily recognized when seen edge-on, the inclination is likely to be near 90 ° and the standoff 
distance, Rq ~ 0.13 pc, similar to bowshocks from stars #2 and #5 in the Kobulnicky et al. 



( 12010( 1 study of Cygnus OB2. 
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3.3. Identification of CO Molecular Complex 

We examined the 12 CO velocity cube o f the IRAS 03063+5735 region from the Canadian 



Galactic Plane Survey (ITaylor et al.ll2003l ). The data cube has velocity channels of width 
0.81 km s _1 and a beamsize of 45' FWHM. There are three peaks of molecular emission 
along this sightline. Figure El shows the 1-dimensional 12 CO spectrum summed over the 
2' square region surrounding IRAS 03063+5735. The first peak near km s^ 1 covers this 
entire region out to 10' from the bowshock nebula with uniform surface brightness. This 
molecular component is likely to be affiliated with very local gas. The second component 
near +5 km s _1 appears as a small filament that crosses near the bowshock but bears no 
morphological resemblance to the diffuse ISM emission seen in the SST and WISE infrared 
images. The third component centered near —40 km s -1 shows a strong channel-by-channel 
morphological resemblance to the mid-IR PAH and hot dust emission. Figure [7] shows a 
color view of the extended IRAS 03063+5735 region with SST [3.6] and [4.5] in blue and 
green, respectively, and the WISE [22] in red. Contours show the 12 CO intensity integrated 

over the —38 42 km s" 1 velocity channels. The 12 CO contours correspond to intensities 

of 5 through 35 K km s _1 in steps of 5 K, equivalent to H2 column densities of 1.5 through 
10.5 xlO 21 cm" 2 ado pting a standard conversion factor of o=3.0xl0 20 cm -2 (K km s -1 ) -1 



(ISolomon et al.lll987l ). Note that the location of the bowshock nebula coincides with a local 
maximum in the H 2 column density of approximately 4.5 xlO 21 cm~ 2 . If we estimate that 
this feature is a roughly spherical molecular cloud with radius r — V = 1.10 pc, then the 
implied number density is Nh2 — 600 cm -3 . We caution that this is only a gross average, as 
12 CO is known to freeze out onto dust grains in the cold dense cloud interiors. 

The similarity of the 12 CO and diffuse ISM morphology, along with the identification of a 
local 12 CO peak at the bowshock location, supports the association of the bowshock with the 
molecular complex. The radial velocity of Star A, which we measure to be —59 ± 18 km s _1 
by a Gaussian fit to the He I A5876 line, adopting a rest wavelength of 5875.6 A. This lends 
further support to the association of the bowshock with this particular molecular complex. 
In principle, one might expect a runaway star to show anomalously large velocities relative to 
the surrounding medium. However, given the clear bowshock morphology of the associated 
nebula, we would also expect such objects to have velocity vectors predominantly in the plane 
of the sky. Hence, the similarity between the radial velocity of Star A and the molecular 
cloud upon which it impinges is consistent with the observed geometry. 



The iBrand fc Blitzl (119931 ) rotation curve in the outer Galaxy shows that radial velocity 
is not a very precise distance indicator along this I =140° sightline. Their Figure 2b shows 
gas at ~ —40 km s -1 associated with a wide range of distances from 2 to 6 kpc. Our 
spectrophotometric distance estimate to Star A of 4.0+1 kpc links the molecular gas at 
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—40 km s 1 to a distance comfortably within that range. 



3.4. IRAS 03063+5735 at other Wavelengths 



IRAS 03063+5735 is not detected to a 3a li mit of 0.037 Jy in the CGPS 1420 MHz 
radio map of this region. IWouterloot et al.l ( 119931 ) searched for but did not detect any OH 
or water masters from this infrared source. 



As first demonstrated by iRubinl (119681 ) the observed flux density of radio emission at 
frequency v (F u ) from an optically-thin, isothermal, H II region at a known distance can be 
related to the H- ionizing photon flux (Q) from the exciting OB star(s). In this study we 
will use values of Q tabulated from model atmospheres of OB stars to determine if an H II 
region surrounding an early B-type star should be visible in our radio continuum images. 



Kerton et al.l ( 119991 ) present a version of th dRubinl (119681 ) formula for observations at 



1420 MHz of an isothermal (75 00 K) H II region containing only hydrogen. Using the stellar 
calibration of ICrowtherl ( 120051 ) we find a B1.5 V star has log(Q) = 46.0 which corresponds 
to F1420 = 900d 2 f, where d is the distance in kpc and / is essentially a covering factor 
(ie., / = 1 would correspond to an ionization bounded region, / = 0.25 could apply to 
a blister- type region). For / = 1 and d = 4 we find F 142 o = 7 mJy This corresponds 
to a source with Tb — 1.1 K in the CGPS ignoring any beam dilution. The average and 
standard deviation of brightness temperature of our 1420 MHz image in a 10x10 arcmin box 
surrounding the position of the bowshock is Tb = 5.68 ±0.08 K so the distance and spectral 
type are consistent with the non-detection of this object in the CGPS radio contiuum images. 
The non-detection in the CGPS image actually puts a fairly strong limit on how early the 
star could be. For example, at 09 V log(Q) = 48.1, and repeating the above calculation 
leads to Tb = 139 K. Even with a very low value for / the H II region associated with a late 
O type star would be easily detectable. All of this is consistant with our non-detection of 
the He II A4511 line in the optical spectrum which shows that Star A is B0. 5 V or later. 



3.5. Bowshock SED 

Figure |8] shows the (background-subtracted) spectral energy distribution of the nebula 
(asterisks) as measured in a large lima-bean shaped aperture in the WISE bandpasses at 
3.6, 4.5, 12 and 22 /mi. The star is detected in the WISE Wl and W2 bands, so we have 
used the higher-resolution SST GLIMPSE photometry to subtract off the star's contribution 
to the WISE Wl and W2 measurements. The diamonds and solid line show the data and 
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extincted blackbody fit for Star A alone. The dashed line shows a iDraine &: Lil (120071 ) dust 
model fit by eye for a radi ation field u = 300 times the average interstellar radiation field 
of 2.17xl0~ 2 erg s _1 cm -2 flMathis et al.lll983T l. a PAH fractional contribution of 1.77% and 
a mean ISM number density of 600 cm -3 . The WISE W2 [4.5] datum lies well above the 
model, possibly because of a contribution from shocked molecular hydrogen in this band. 
The IRAS 60 and 100 /xm data also show relatively poor agreement with the model, but 
the uncertainties on these data are large owing to the large (4') beamsize and background 
subtraction complications. The diffuse infrared surface brightness varies substantially over 
the nebula, suggesting substantial spatial variations in the local dust and gas density. Vari- 
ations in temperature and density throughout the nebula are certain to render the single 
averaged SED a gross oversimplification. Hence, the dotted line fit is mostly for illustrative 
purposes, especially given large uncertainties on the local number density in the vicinity of 
the bowshock. 



3.6. Search for Associated Stellar Cluster 

The apparent increase in stellar density just to the north-east of the Star A initially 
suggested the possibility of a stellar cluster at this location. However, a color-magnitude 
diagram using JHK photometry of 2MASS Catalog sources within a 1' radius shows that 
the stars have a considerable spread in color (hence, reddening), and no clustering in either 
color-color or color-magnitude diagrams is evident. Therefore, we find no evidence of a stellar 
cluster in this vicinity, at least to the depths probed by the 2MASS and GLIMPSE catalog 
data. 



4. Discussion 

4.1. Runaway or Low- Velocity Star? 

The expected proper motion of an object at 4.0 kpc distance with a transverse motion 
of 30 km s -1 is ~1 mas yr _1 , too small to have a measured proper motion even over decade 
time scales in multi-epoch photographic plate surveys. Indeed, the USNO-B1.0 Catalog does 
not report a proper motion. Hence, there is no direct evidence for Star A having a high space 
velocity. Rather, it is inferred as a possibility from the morphology of the nebula. 



It is also possible that low-velocity stars moving in a medium with a density gradient 
(jWilkinl 120001 ) can produce a bowshock-like nebula. The proximity of Star A to a molecular 
cloud having a maximum column density just ahead of the bowshock (see Figure [7]) makes 
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this explanation worth considering. Perhaps Star A formed within or near the pictured 
molecular complex. However, formation of such a massive star typically requires relatively 
high column densities characteristic of massive molecular clouds and would be accompanied 
by the formation of other lower mass stars comprising a small embedded cluster. We find no 
evidence of either high column density molecular clouds (the highest is 10.5 x 10 21 cm" 2 at 85" 
=1.6 pc north of the bowshock) or a stellar cluster containing lower mass stars. Furthermore, 
the observed excess emission in the [4.5] band (Figure [8]) would be well explained by emission 
from shocked H2 if the nebula were a genuine shock. It is not clear whether a low-velocity 
B1.5 star carving a cavity within a local density gradient would generate a wind shock 
sufficient to produce excitation of H 2 . More detailed spectral information and shock modeling 
would be needed to determine if this is a possibility. Given the current data, we prefer the 
interpretation that Star A is encountering the molecular cloud for the first time after having 
traveled some distance from its point of formation at a high velocity. 



4.2. Possible Origins of a Runaway Star 

The orientation of the cometary bowshock nebula suggests that Star A's direction of 
travel is roughly toward the equatorial north. We searched the SIMBAD database for young 
stars, star clusters, and H II regions up to 3 degrees to the south to identify candidate 
birthplaces of Star A. At a distance of 4.0 kpc and an assumed tangential velocity of 30 
km s -1 (30 pc Myr -1 ) in the plane of the sky, it would be reasonable to suppose that the 
star has traversed up to 600 pc (8°on the sky ) or even more during its fewxlO Myr main- 
sequence lifetime. However, given that the star is already below the Plane and moving 
toward Galactic north, adopting the maximum travel time allowed by the main-sequence 
lifetime implies a formation locale well outside the scale height of Milky Way molecular 
clouds required to form such a massive star. Therefore, we prefer the interpretation that 
the runaway has been ejected recently and confine our search within several degrees of its 
current location. 

One likely candidate for Star A's birthplace is the H II region and stellar cluster associ- 
ated with IRAS 03064+5638 located 0.94° (64 pc) due south. Figure [9] shows a three-color 
image of the field around Star A with the CGPS 1420 MHz radio continuum in red, the WISE 
[11.1] in green, and the WISE [4.5] in blue. Contours indicated the 12 CO brightness of 5, 10, 
15, 25, 30 K km s -1 integrated between LSR velocities of —31 and —43 km s -1 . The bright 
ridge of radio emission extending 2° north to south is a thermal rad io filament at an un- 



known distance and powered by a yet-unidentified source (lGreenlll989l ). Both the bowshock 
nebula IRAS 03063+5735 and IRAS 03064+5638 have 12 CO peaks in this velocity range, 
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sug gesting a comm o n dist ance. IRAS 03064+5638 contains a stellar cluster characterized 



by ICarpenter et al.l (119931 ) as having a few mid-B type stellar members, for their adopted 
distance of 2.2 kpc. If located at 4.0 kpc instead, these become early B stars of comparable 
mass and age to Star A. The 1420 MHz radio flux of IRAS 03064+5638 is 101 ± 3 mJy, 
equivalent to a B0 V star at the adopted distance. Apparently, the IRAS 03064+5638 clus- 
ter contains at least one star more massive than Star A, making it an environment from 
which Star A could plausibly have been ejected. Three-body encounters which eject one star 
often result in the formation of a tight binary. Therefore, our working scenario predicts the 
existence of close early-B binary within IRAS 03064+5638. However, the extinction appears 
high at this location, and no suitably bright optical sources are visible. 

The peak molecular column density at IRAS 03064+5638 is approximately A r (if 2 )=4.5xl0 21 
cm -2 . The adjoining molecular clouds which lack bright IR emission have peak N(H 2 ) al- 
most twice this value, suggesting other possible sites of intermediate-mass star formation 
nearby. These, yet unidentified, sites of star formation could have also been the birthplace 
of star A. 

Another possible origin for Star A is the open cluster NGC 1220, located a full 3° (~201 



pc) to the south of Star A, but along the same vector pictured in Figure [9j lOrtolani et al. 



(120021 ) identified 26 probable members of this open cluster, estimating a distance of 1800+200 
pc and an age of 60 Myr. While this denser and more populous cluster provides a suitable en- 
vironment for ejecting a B star via N-body encounters, the much smaller heliocentric distance 
(if correct) makes a physical connection between NGC 1220 and IRAS 03063+5735/Star A 
less likely. 



5. Conclusions 



We have identified an IRAS infrared source IRAS 03063+5735, previously reported as 
a galaxy, to be a bowshock nebula powered an early B-type (B1.5 V) star. The infrared 
source was selected from among several hundred sources having IRAS color-color criteria 
matching those expected from soft UV radiation fields generated by intermediate- mass (i.e., 
B-type) stars. Our measurement of the —58 km s _1 LSR radial velocity of the star, combined 
with 12 CO line maps, allow us to identify a molecular cloud complex at a common velocity 
and ~4.0 kpc distance. The bowshock is plausibly produced by the supersonic motion of a 
runaway B star as it impinges upon a local maximum in the molecular column density of 
~ 4.5 x 10 21 cm -2 . However, a low-elocity star within a local density gradient may also be 
capable of producing the observed phenomena. We identify the H II region and stellar cluster 
associated with IRAS 03064+5638 at a projected distance of 64 pc as one plausible birth 
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site for a runaway star. The association of a star having a spectrophotometric distance with 
a particular molecular feature adds another data point to help anchor the (poorly defined) 
rotation curve in the outer Galaxy. 
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Table I. Star A and Comparison Star Line Measurements 



Star EW(HeII A5411) EW(HeI A5876) EW(Ha) Spectral Type Ref. 

(A) (A) (A) 



Star A 


0.11 (0.10) 


0.78 


(0.11) 


3.70 


(0.71) 


B1.5 V 


1 


DL Cam 


0.24 (0.01) 


1.03 


(0.06) 


2.20 


(0.34) 


BO III 


2 


HD36822 


0.13 (0.01) 


0.84 


(0.01) 


2.52 


(0.31) 


BO III 


2 


HD23180 


<0.01 


0.78 


(0.03) 


2.34 


(0.43) 


Bl III 


2 


HD30836 


<0.01 


0.80 


(0.02) 


2.17 


(0.24) 


B2 III 


2 


HD29248 


<0.01 


0.78 


(0.02) 


2.44 


(0.31) 


B2 III 


2 


HD31237 


<0.01 


0.81 


(0.02) 


2.28 


(0.27) 


B3 III 


2 


e Cas 


<0.01 


0.52 


(0.01) 


3.37 


(0.52) 


B3 III 


2 


S Per 


<0.01 


0.42 


(0.01) 


3.87 


(0.51) 


B5 III 


2 


A Cet 


<0.01 


0.32 


(0.01) 


4.69 


(0.68) 


B6 III 


2 


HD16727 


<0.01 


0.22 


(0.01) 


4.79 


(0.74) 


B7 III 


2 


7 Per 


<0.01 


0.19 


(0.01) 


4.60 


(0.65) 


B8 III 


2 


MT378 


<0.09 


0.81 


(0.05) 


2.19 


(0.14) 


BO V 


3 


MT429 


0.13 (0.03) 


0.69 


(0.02) 


2.43 


(0.18) 


BO V 


3 


HD22951 


0.02 (0.01) 


0.71 


(0.02) 


2.59 


(0.32) 


B0.5 V 


2 


MT605 


0.05 (0.02) 


0.64 


(0.02) 


2.31 


(0.15) 


Bl V 


3 


MT365 


0.05 (0.02) 


0.72 


(0.03) 


2.43 


(0.27) 


Bl V 


3 


MT605 


0.05 (0.02 


0.64 


0.02 

V / 


2.31 


(0.15) 


Bl V 


3 


MT187 


0.03 (0.02) 


0.67 


(0.02) 


2.39 


(0.16) 


Bl V 


3 


HD21856 


<0.01 


0.68 


(0.02) 


2.61 


(0.33) 


Bl V 


2 


HD24131 


< 0.01 


0.66 


(0.01) 


2.74 


(0.36) 


Bl V 


2 


HD11241 


<0.01 


0.81 


(0.02) 


2.95 


(0.56) 


B1.5 V 


2 


MT234 


<0.01 


0.67 


(0.03) 


2.90 


(0.27) 


B2 V 


3 


HD23990 


<0.01 


0.67 


(0.01) 


2.11 


(0.41) 


B2 V 


2 


MT298 


<0.01 


0.64 


(0.07) 


3.44 


(0.54) 


B3 V 


3 


29 Per 


<0.01 


0.60 


(0.01) 


3.82 


(0.56) 


B3 V 


2 


34 Per 


<0.01 


0.55 


(0.01) 


4.04 


(0.63) 


B3 V 


2 


53 Per 


<0.01 


0.58 


(0.01) 


3.83 


(0.54) 


B4 IV 


2 


31 Per 


<0.01 


0.58 


(0.02) 


3.77 


(0.65) 


B5 V 


2 


HD21278 


<0.01 


0.47 


(0.01) 


4.23 


(0.63) 


B5 V 


2 
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Table 1 — Continued 



Star 


EW(HeII A5411) 


EW(HeI A5876) 


EW(Ha) 


Spectral Type 


Ref. 




(A) 


(A) 


(A) 






HD24504 


<0.01 


0.43 (0.02) 


4.31 (0.76) 


B6 V 


2 


HD18537 


<0.01 


0.23 (0.01) 


4.84 (0.80) 


B7 V 


2 


P Per 


<0.01 


0.23 (0.01) 


4.30 (0.65) 


B8 V 


2 



Note. - 1-This work; 2-SAO Star Catalog; 3- lMassev fc Thompson! fll99lh : 
Kiminki et all f|2007h 
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Table 2. 2MASS 03101044+5747035 Photometry 



Wavelength 


Mag. 


Flux 


Source 


(fim) 




(mJy) 




0.44 


16.33 


1.25 


1 


0.65 


14.54 


4.09 


1 


0.75 


13.93 


5.68 


1 


1.2 


12.04 


26.2 


2 


1.6 


11.38 


28.4 


2 


2.2 


11.03 


24.0 


2 


3.6 


10.82 


13.8 


3 


4.5 


10.78 


9.19 


3 



Note. - 1-The USN O-Bl.O Cat- 
alog. JMonet et all l2003h: 2-2M ASS 
(Th e Two Micron All Sky Survey 



20061 ): 3-GLIMPSE 360 Point Source 
Catalog 
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Table 3. IRAS 03063+5735 Photometry 



Wavelength 


Flux 


Notes 


(/xm) 


(Jy) 




Wl 3.6 


0.21 


1 


W2 4.5 


0.23 


1 


W3 12.1 


4.6 


1 


W4 22.0 


4.9 


1 


IRAS 60 


33.7 


2 


IRAS100 


83.8 


2 



Note. - - We adopt a 10% 
uncertainty on the flux mea- 
surements at all WISE band- 
passes and a factor of 2 for 
the two IRAS measurements. 
1-The reported flux reflects 
the entire nebula in a large 
aperture minus the flux from 
the star, as measured from 
SST GLIMPSE photometry, 
as reported in Table [2j 2- 
IRAS Point Source Cata- 
log value. The IRAS mea- 
surements at 12 and 25 
are systematically a factor 
of about 1.8 lower than the 
W3,W4 WISE measurements 
in similar bandpasses. This 
discrepancy probably reflects 
the much larger IRAS beam- 
size and associated uncer- 
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tainty in subtracting the dif- 
fuse background emission. 
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Fig. 1.— (right panel) SST three-color image in Galactic coordinates of the 
IRAS 03063+5735 region with SST [3.6] in blue, [4.5] in green, and WISE W3 12 fjm in 
red, revealing the double bowshock morphology of the nebula, (left panel) 2MASS J/H/K 
image of the same field in blue/green/red, respectively. The star 2MASS03101044+5747035 
(Star A) lies at the center of the nebula and can be seen in both images. 




Fig. 2.— Averaged WIRO spectrum of 2MASS03101044+5747035. Locations of stellar He II 
A5411 (not detected), He I A5876 and Ha are marked. Other absorption lines are interstellar, 
including Na I AA5889,5995 and several diffuse interstellar band features. 
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Fig. 3. — Ratio of He II A5411 equivalent width to He I A586 76 EW versus spectral 
class/temperature for early ty pe stars from the [Jacoby et al.l (11984 J84) spectral atlas (up- 
per ordinate) and the Tlusty ( Lanz fc Hubeny 20031 ) model atmospheres (lower ordinate), 
for both luminosity class V (dwarf) and III (giant) stars, as indicated by the legend. The 
dashed line shows the third-order polynomial fit to the luminosity class V models and spec- 
tral atlas stars, as described by the equation. Star A has a ratio EW(HeII A5411)/EW(HeI 
A5876)< 0.14, indicating a spectral type later than BO. 
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Fig. 4. — (upper panel) EW of He I A5876 versus temperature/spectral type for B type stars. 
(lower panel) Ratio of Ha to He I A5876 equivalent width versus temperature/spectral type. 
Pluses denotes dwarfs (Lum. V) fro m Table ED while asterisk s denote giants (Lum. III). 
Triangles are model atmospheres from lLanz fc Hubenyl ( 120071 ) . The dashed curve shows a 
3rd-order polynomial fit to the pluses and asterisks. The solid line and dotted lines show 
the measurements and la uncertainties for Star A. 
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Fig. 5. — Spectral energy distribution for Star A. The curve is a blackbody with T=23,000K 
and radius R=6.1 a ppropriate for a B1.5 V star, interpolated f rom the data of 
Drilling fc Landoltl (2000), with Ay = 5 .4 mag and a Cardeili et al. ( 1989 ) extinction law in 



the optical and llndebetouw et al.l (120071 ) extinction law in the mid-IR. A distance of 4.0 kpc 
yields a good fit. 
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Fig. 6.— A 1-dimensional 12 C0 spectrum from the CGPS along the IRAS 03063+5735 
sightline averaged over a square 2' region. The molecular component near -40 km s _1 bears 
a strong morphological similarity to the diffuse emission from warm dust. 
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Fig. 7. — A color representation of the larger IRAS 03063+5735 region in Galactic coordi- 
nates with SST GLIMPSE [3.6] in blue, [4.5] in green, and the WISE [22] in red. Contours 
show the 12 CO brightness map from the CGPS integrated over the LSR velocity range -38 
to -42 km s _1 . The 12 CO contours correspond to column densities of 1.5 xlO 21 cm~ 2 - 
10.5 xlO 21 cm -2 in steps of 1.5 xlO 21 cm -2 . 
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Fig. 8. — The spectral energy distribution of the IRAS 03063+5735 bowshock nebula (aster- 
isks) from the four WISE bandpasses and two IRAS measurements at 60 and 1 00 /im. Dia- 
mond s and the solid line show the extincted stellar SED. The dotted line shows a lDraine fc Li 
(200 7J) dust model w ith radiation intensity u = 300 times the mean interstellar radiation field 
jMathis et al1ll983l ) and PAH contribution 1.77%. 
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Fig. 9. — Color representation of the field around IRAS 03063+5735 as seen in the 1420 
MHz CGPS radio continuum map (red), WISE W3 [11.1] (green), and WISE W2 [4.5] 
(blue). Contours depict the 12 CO integrated between LSR velocities of -31 km s -1 and 
-43 km s _1 . A green arrow points rearward along the bowshock axis of symmetry to a 
hypothetical birthplace of Star A near the H II region and stellar cluster associated with 
IRAS 03064+5638 at a projected distance of 65 pc. 



